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ABSTRACT 

Following our previous work, we investigate through hydrodynamic simulations the destruction 
of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an 
idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of 
Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. 
We vary the relative velocity between the reverse shock and eject a clump to explore the effects of 
shock- heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider 
gas metallicities from Z/Z e = 1 to 100 and their influence on cooling properties of the cloud and the 
thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate 
grain sputtering for a variety of species and size distributions. In the metallicity regime considered 
in this paper, the balance between increased radiative cooling and increased grain erosion depends on 
the impact velocity of the reverse shock. For slow shocks (v s hock < 3000 km s _1 ), the amount of dust 
destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. 
For higher shock velocities (v s hock > 5000 km s _1 ), an increase in metallicity from Z/Zq = 10 to 100 
can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction 
varies widely across grain species and simulation parameters, our most extreme cases result in complete 
destruction for some grain species and only 44% dust mass survival for the most robust species. These 
survival rates are important in understanding how early supernovae contribute to the observed dust 
masses in high-redshift galaxies. 

Subject headings: hydrodynamics — supernova remnants — shock waves — dust 



1. INTRODUCTION 

Over the last two decades, far-infrared (FIR) and mil- 
limeter observations of high-redshift quasars (z > 6) have 
produced estimates for gal actic dust masses a s high as 
1Q 8 Ma) ([ S mail et al.|[l997} | Hughes etldll998| |Bertoldi 
|et al.||2003[ j Wang et al^^p . In order to explain the 
formation of this large quantity of dust within the short 
lifetime (~1 Gyr) of the universe at this epoch, a mecha- 
nism must exist which is capable of both significant and 
rapid dust production. One recently pursued solution 
is that the majority, if not all, of this dust comes from 
core-collapse supern ova (CCSN) explosions of the first 
enerations of stars ([M organ fc Edmunds|2003 Maiolino 



et al. 2004; Hirashita et al. 
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ch neider j |2007| jNozawa et al.j 2010) has suggested 
0.1-0.3 M e of dust could be formed per super- 



1.0 Mq required to explain high-redshift dust (Morgan 
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Dwek et al.||2007| |Dwek & (Jherchnetl 
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2011). However, observational efforts 



orders of magnitude short of these values ( Stanimirovic 
et al.||2005| [Williams eFaLl[2006l IMeikle et al.||2007MRh( 
eT"aLl[2008l l2009HKotak et al ||2009| ). More recent obser 
vations aimed at finding colder dust (T < 40 K) using 
AKARI and the Balloon-borne Large Apera ture Submil- 



limeter Telescope (B LAST) by|Sibthorpe et ah 
well as Herschel by |Barlow et al. ( |2010[ ) aricT 



(201(| 
Matsuura 



as 



et al. (2011), find larger dust masses. Matsuura et al. 
Wn\ 15^0.6 M© of dust in SN1987A. Such estimates 
could make the argument for SNe as dust factories more 
plausible. 

One question that arises from these SNR studies is 
what fraction of the freshly formed dust predicted by 
theory will survive the interaction between the reverse 
shocks and the ejecta. As the fast-moving dust (V e j > 
1000 km s _1 ) is impacted by the reverse shock, it will 
be subject to sputtering and grain-grain collisions as the 
density and temperature of the dust-enriched gas are in- 
creased. 



In our previous work (Silv ia et al.||2010 hereafter Pa- 
per I), we investigated, through hydrodynamic simula- 
tions, the destruction of newly formed dust grains by 



crushing" simulations (Woodward 


1976; 


1994 


Klein et al.|1994 


) , we found that tl 
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bution of the dust grains as well as the initial density 
of the ejecta cloud and relative velocity between the re- 
verse shock and the cloud. In the most extreme cases, we 
found grain destruction to vary from 20-100% depending 
on grain species. We also found morphological similari- 
ties within these simulations to the observationa l studies 
of Cassiopeia A presented in |Fesen et al.| ( |2011[ ) , specif- 
ically in the scenarios where the effects of cooling led to 
the fragmentation of ejecta clouds. 

However, our earlier simulations computed the radia- 
tive cooling and sputtering rates for approximately so- 
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lar metallicity (Z@) gas. In reality, the gas contained 
in supernova ejecta will be highly metal-enriched. The 
question arises as to the balance between the enhanced 
radiative cooling and the increase in sputtering yields in 
this high-metallicity regime. This balance will also be 
influenced by the changing density of the metal-enriched 
gas as the cloud is shredded and the plasma thins out. 

In the current paper, we perform additional cloud- 
crushing simulations to study the evolution of dust mass 
contained within an idealized ejecta knot as it is im- 
pacted by a supernova reverse shock. In Section [2] we 
give a brief review of the code used to carry out these 
simulations, the methods for tracking our dust popula- 
tions, and the changes made from our previous work. In 
Section [3j we describe the simulations unique to this pa- 
per, specifically aimed at probing the higher metallicity 
regimes expected to be present in supernova ejecta. We 
present the results of these simulations in Section [4] and 
conclude in Section [5] with a summary and discussion of 
the implications of these results. 

2. METHODOLOGY 

2.1. Code and Simulation Setup 

As in Paper I, we use the Eulerian adaptive mesh re- 
finement (AMR), hydrodynamics + N-body code, Enzo 
(|Bryan fc Norman||1997||Norman fc Bryan||1999| |0'Shea 



et al.|2UU5| ). We cfo not make use ol any ol the cosmolog 
ical or gravity- solving components of the code owing to 
the idealized nature of our problem. We follow the same 
cloud-crushing setup as outlined in Paper I and refer the 
reader to that work for a detailed description. 

The user-supplied parameters required to initialize the 
simulation are: cloud radius, r c i ou d; cloud temperature, 
^cioud? density of the ambient medium, p m ; initial over- 
density, °f the cloud with respect to the ambient 
medium (p c ioud = XPm); and the velocity of the shock 
relative to the stationary cloud, v s hock- From these input 
parameters, the following values must be derived in order 
to completely initialize the simulation: temperature of 
the ambient medium, T m , post-shock density, p s hock> and 
post-shock temperature, T s h oc k- We set T m so that the 
cloud remains in pressure equilibrium, while the shock- 
related values are calculated using the Rankine-Hugoniot 
jump conditions. 

During runtime, we take advantage of the AMR capa- 
bilities of Enzo by employing the same refinement cri- 
terion as described in Paper I, increasing resolution in 
areas of the simulation that contain a significant fraction 
of cloud material. Cells that are initially enclosed within 
the cloud radius are assigned a "cloud material" value, £, 
that is advected with the flow in the same way as density. 
When a cell exceeds a pre-defined cloud material "mass" , 
= £V^eiij where T^eii is the cell volume, the resolution 
of that cell is doubled. A more thorough description of 
this refinement process can be found in Paper I. 

2.2. Dust Tracking and Post-processing 

Our dust-grain populations continue to be tracked 
through tracer particles embedded in the flowing gas. 
As before, we post-process the density and temperature 
histories of each tracer particle to compute dust survival 
rates. For the initial distributions in grain radii, we again 
made use of the values calculated by|Nozawa et al. (2003) 



for a CCSN with a progenitor mass of 20 M (see Fig- 
ure [TJ. We follow the evolution of all nine grai n species 
includ ed in the unmixed grain model of |Nozawa et al. 
(2003). To determine the evolution in dust mass, we use 
the same mass proxy outlined in Paper I, which involves 
tracking the changes in n(a) x a 3 , where n(a) is the num- 
ber density of grains and a is the grain radius. However, 
the erosion rates used to sputter the grains with each suc- 
cessive time step are different than those used in Paper I, 
as described in the following section. We note that we 
only account for dust grain destruction in this work and 
do not include a prescription for possible grain growth 
that might occur at high densities and low temperatures. 
For the majority of our simulations, the amount of time 
that the dust grains spend in an environment conducive 
to grain growth is relatively brief and any potential in- 
crease in dust mass should be minimal. 




a [cm] 

Fig. 1. — Reproduced grain radius distributions for the species 
expected in the unmixed ejecta model of a core-collapse supernova 
with a progenitor mass of 20 Mq as calculated by |Nozawa et aT] 
(2003). The y-axis is the abundance of grains of a given radius, 
J~(a) : such that the number density of grains between a and a + 5a 
is n(a) = f(a)Sa, where 5a is set by the number of bins used to 
track the distribution. 



2.3. Modifications and Additions 

To account for the fact that we expect the super- 
nova ejecta clumps to be metal-enriched, we implement 
both new cooling curves and grain erosion rates for high- 
met allicity__gas : __WecaIciilate the cooling rates using 
Cloudy ( |Ferland et al. 1998[ ) for gas assumed to be in 
ionization equilibrium. In addition, these rates assume 
that both the electrons and the different ion species have 
the same temperature. Cooling due to the thermal emis- 
sion of dust grains is not included. 

For this work, we compute rates for metallicities of 
Z/Zq = 1, 10, and 100. The element abundances for 
solar metallicity gas come from the composition listed 
in the documentation for Cloudy, where all abundances 
are specified by number relative to hydrogen. We de- 
fine the metallicities with values greater than unity to 
mean that the abundances for metals are increased by 
factors of 10 and 100 from their solar values. The abun- 
dance of hydrogen and helium remain the same for all 
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three metallicities. To simplify the notation, we refer to 
these metallicities as 1 Z , 10 Zq, and 100 Zq for the 
remainder of the paper. We also note that the values 
for the metal abundances remain static for the duration 
of our simulations, and any potential increase in metal 
abundance due to sputtered dust grains is omitted. We 
comment on this omission in the final section of the pa- 
per. Figure [2^i shows the cooling rate coefficients for 1, 
10, and 100 Zq gas as a function of temperature. To un- 
derstand the contributions that hydrogen, helium, and 
heavier elem ents make on the cooling curves, we refer 
the reader to |Gnat fc Sternberg! ( |2007[ ). 

In contrast to Paper 1, in which the cooling rates were 
applied to all cells within the simulation domain, we only 
cool those cells that contain cloud material. This reduces 
the computational time required to calculate the cooling 
rates and accounts for the fact that, while the cloud itself 
can be highly metal-enriched, the ambient medium is at 
much lower densities and metallicities. To ensure that 
important information is not lost by cooling only the 
cloud, we consider that for an ambient medium with a 
metallicity of Z = 1 Zq, a number density of n = 1 cm -3 , 
and a temperature of T = 10 6 K, the cooling time is a 
few x 10 4 yrs, much longer than the time scale for any 
of the simulations included in this work. 

We ca lculate new erosion r ates using the formula pro- 
vided by Nozawa et al. (2006) for the limit in which ther- 
mal sputtering is the dominant mechanism for grain ero- 
sion. This is appropriate, since the tracer particles used 
to track the dust are embedded in the flow. In com- 
puting these erosion rates, we use the same scaled solar 
metal abundance ratios that were used to produce our 
new cooling functions. Computing rates based on these 
enhanced metal abundances is an important step. The 
sputtering yield, Y(E) where E is the energy of the im- 
pacting ion, is strongly dependent on the atomic mass 
of the impactor; high-mass ions can have orders of mag- 
nitude higher yields at high energies. The differences in 
erosion rates as a function of metallicity are shown for 
carbonaceous, silicate, and ferrous grains in Figure [2}d. 
The differences do not become significant until the tem- 
perature exceeds a few times 10 6 K and, even then, only 
for the highest metallicity. In order to simplify Figure 
[2}d, we have omitted the other six grain species studied 
in this work. However, the erosion rates for the omitted 
species are comparable to the three presented species. 

We include Figure (3^i to show the nature of the sput- 
tering yields for a sample of impacting ions as a function 
of energy and Figure [SJd to show the contribution of vari- 
ous impacting ions to overall grain sputtering for gas with 
both a metallicity of 1 Zq and 100 Zq. For simplicity, we 
only present these quantities for C grains, but note that 
the figures are qualitatively similar for the other grain 
species. From Figure (3Jd, it is evident that the contribu- 
tions from the sputtering yields of high-mass ions do not 
begin to dominate over hydrogen and helium until the 
metallicity reaches 100 Zq, which agrees with the trend 
observed in the erosion rates presented in Figure^. 

When these new cooling and sputtering rates are im- 
plemented within our simulations, we find that, for the 
high metallicity (Z = 100 Zq) simulations, the densest 
cells end up with extremely short cooling times, often 
shorter than the time step set by the Courant-Friedrichs- 
Lewy (CFL) condition. This tends to result in cells that 
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(a) Metallicity-dependent cooling curves 
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(b) Metallicity-dependent erosion rates 

(a) Radiative cooling curves for solar-scaled metal abun- 
= 1 Zq, 10 Zq, and 100 Zq. (b) Dust grain erosion rates 



for C, SiC>2, and Fe grains. 

over-cool during a given time step and can lead to nega- 
tive cell energies. To avoid this issue, we make two mod- 
ifications to Enzo. First, we add a temperature floor, 
Tfi oor = 1000 K, to all of our simulations to prevent 
clouds from becoming unrealistically cold. Given the en- 
ergetic environment of the supernova remnant and the 
background radiation from shock-heated, x-ray emitting 
gas in the remnant's shell, such a floor is physically rea- 
sonable. Second, we modify the time-step calculation 
such that it is never more than 25% of the cooling time. 
While this prevents negative cell energies, it can become 
computationally expensive, as the cloud is compressed to 
high densities and the cooling time becomes very short 
in some regions of the cloud. It is this computational 
limitation that prevents us from exploring even higher 
metallicities, at least within the bounds of our computa- 
tional resources. 

3. SIMULATIONS 

All of the simulations presented in this work have root- 
grid dimensions of 512 x 256 x 256 and a physical res- 
olution of 1.25 x 10 15 cm (~84 AU) per cell edge. As 
in our previous work, we initialize each cloud to have a 
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TABLE 1 

Simulation Parameters 



i(r i(r 
E (eV) 

(a) Sputtering yields 
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E (eV) 

(b) Metallicity-dependent contributions to sputtering 



10* 



Fig. 3. — (a) Sputtering yields, Y{(E), of C grains versus pro- 
jectile energy, E, for various incident ions; (b) The contribution of 
various incident ions to the sputtering of dust grains, Qi, where Qi 
is defined to be the product of the element abundance by number 
relative to hydrogen, ni/riH, and the sputtering yield normalized 
to that of hydrogen, Y[/Yn, as a function of energy. We show 
the contributions based on solar abundance ratios (solid lines) and 
abundance ratios for Z = 100 Zq (dashed lines). 

radius of r c i ou d = 10 16 cm, which spans 8 cells at the 
root-grid resolution. We also use the same Gaussian en- 
velope formulation as in Paper I, with the density fall-off 
occurring at r = 0.7 r c i ou d. 

We allow the simulation domain to be refined up to 
three additional levels such that the highest resolution 
cells will be 8 times smaller than the root-grid. This dif- 
fers from our previous work, where we only allowed for 
two additional levels of resolution. The third level of re- 
finement was required specifically in the simulations with 
high metallicity, to follow the collapse of some fragments 
to much higher densities and smaller spatial scales than 
in the low- metallicity cases. When we allow for a fourth 
level of refinement, the final dust masses do not change 
significantly; we therefore limit ourselves to three levels 
to save computational resources. As mentioned above, 
our cell refinement is based on the amount of cloud ma- 
terial in a given cell. 

For this work, we primarily focus our exploration of 
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Note. — Parameter definitions: % is the initial over- 
density of the cloud with respect to the ambient medium, 
v s hock is the relative velocity between the inflowing shock 
and the stationary cloud, and t cc is the cloud-crushing time. 

parameter space to the relative velocity between the re- 
verse shock and the eject a cloud and the metallicity of 
the gas contained within the cloud. Specifically, we in- 
vestigate shock velocities of 10 3 km s _1 , 3xl0 3 km s _1 , 
5xl0 3 km s _1 , and 10 4 km s _1 and metallicities of 1 Z , 
10 Z , and 100 Z@. For ease of reference, we provide 
a numbered list of all simulations in Table [TJ For all 
simulations, the cloud constrast is x = 1000, the initial 
cloud temperature is set equal to our temperature floor, 
^cioud = ^floor = 1000 K, and the total runtime is 4.2£ cc , 
where t cc = X^Tcioud/^shock is the cloud-crushing time 



flKlein et al.|1994| ). 

Figure [3] snows the evolution in density of the shock- 
cloud interaction from the initial impact to t ~ 3 t cc 
for v s hock = 5 x 10 3 km s _1 and Z = 100 Z Q . Notable 
features include the tail of ablated material produced as 
the shock washes over the cloud, and the numerous cold, 
dense fragments that form a result of the high cooling 
rates for the metal-enriched gas. A handful of these frag- 
ments persist for a large fraction of the total simulation 
time. 

While our previous work ran all simulations for a total 
of ten cloud-crushing times, we found that in our new 
Simulations 2 and 3 (slow shock and high metallicities) 
it was difficult to keep the majority of the dust parti- 
cles within the computational domain much beyond four 
cloud-crushing times without significantly increasing the 
size of the root grid. Though possible, we deemed this so- 
lution to be a poor use of our computational resources. 
As will be seen in Section [4j the dust masses are still 
evolving in the relatively slow-shock (v s hock = 10 3 km s _1 
and 3xl0 3 km s _1 ) simulations at the point of ter- 
mination. However, the fast-shock simulations do not 
have significant changes in dust mass beyond that time. 
Therefore, while we do see clear trends in the slow-shock 
simulations as a function of metallicity, we cannot say 
definitively how these trends would behave at later times. 

4. RESULTS 

We present the results for the twelve simulations listed 
in Table ffl with a focus on the evolution in total dust 
mass and the time spent by the post-processed dust 
particles in various areas of density-temperature phase- 



Reverse Shock Dust Destruction 



-23.2 



-22.4 



logio[Density (g cm -3 )] 
-21.6 



-20.8 



1 1 1 1 






0.41 


0.81 


W% ■ 

1.22 


mm 

1.62 


2.03 


2.43 


2.84 


3.24 



4.1. Dust Mass Evolution 



Fig. 4. — Density projections in the x-z plane weighted by cloud material for the case of ^shock = 5000 km s _1 and Z = 100 Zq. The 
number in the lower left corner of each panel indicates the ti me of the sna pshot in units of the cloud-crushing time, t cc - Projections were 
made using the software analysis package, yt (yt-project. org; [Turk et al.|2011[ ). 

space. crease in dust destruction (see the results of Simula- 

tion 3). In the case of v s hock = 3 x 10 3 km s _1 and 
Z = 100 Zq, the amount of dust destruction is noticeably 
less for some grain species (AI2O3, MgSiOs) and consid- 
erably more for others (Fe, Si). This suggests that, for 
some grain species, the increased cooling due to the en- 
hanced metallicity is able to counter-balance the shock- 
heating and keep the gas in a regime of decreased sput- 
tering. For the grain species with lower survival rates, a 
shock velocity of 3 x 10 3 km s _1 is sufficient to drive up 
the amount of dust destruction. However, we again note 
that this is only based on the results at t = 4.2£ cc . At this 
velocity, the simulations with Z = 1 Zq and 10 Zq ap- 
pear to have hit a dust mass plateau, but for Z = 100 Zq 
the evolution in dust mass still has a clear downward 
slope. For v s h oc k = 5 x 10 3 km s _1 and 10 4 km s _1 , 
an increase in metallicity leads directly to an increase in 
total dust destruction. The most drastic change occurs 
when the metallicity is increased from Z = 10 Zq to 
100 Z Q ; ^20% more of the initial dust is lost with this 
change in metallicity for Fe and Si grains in the high- 
est velocity scenario. For all simulations in which the 
evolution in dust mass plateaus at late times, the cloud 
has been significantly shredded and the density of the 
gas has been reduced to a level that does not allow for 
appreciable grain sputtering. 

In studying the final (surviving) dust masses for each 
grain species, we find that the values are widely scat- 
tered for the simulations with the slower shock velocities. 
The higher velocity simulations show dust survival of less 
than 20% for six of the nine species, with C, Fe, and Si 
grains being the three species with the least destruction. 
Although Fe grains had a higher survival rate than other 
dust species, the amount of Fe dust destruction was con- 
siderably higher than any of the simulations in our pre- 



The evolution of dust mass for each of the nine dust 
species and for all twelve simulations can be seen in Fig- 
ure [5] and the final dust masses can be located in Ta- 
ble [3 The most notable difference between the sim- 
ulations is the drastic increase in dust destruction be- 
tween the slowest shock velocity, v s hock = 10 3 km s" 1 , 
and the highest shock velocity, v s hock = 10 4 km s _1 , 
in some cases a difference of >70%. For the interme- 
diate shock velocities, v s hock = 3 x 10 3 km s _1 and 
5 x 10 3 km s _1 , the evolution in dust mass depends on 
the grain species. For many species, the simulations with 
^shock = 5 x 10 3 km s _1 show comparable dust destruc- 
tion to the highest velocity shocks, while the simulations 
with v s hock = 3 x 10 3 km s _1 exhibit only moderate dust 
destruction. There does appear to be a fundamental dif- 
ference between the simulations with the slowest shock 
velocities and those with faster shocks at early times. 
For nearly all grain species, there is a sharp drop in dust 
mass within one cloud-crushing time for shock velocities 
of 3 x 10 3 km s _1 and higher, with the magnitude of the 
drop increasing with shock velocity. This suggests that 
there is a threshold velocity somewhere between 10 3 and 
3 x 10 3 km s _1 at which the initial impact of the shock 
into the cloud results in a degree of shock heating and 
compression that significantly influences the overall dust 
destruction in the cloud. For slower shocks, the majority 
of the dust destruction occurs at much later times, as 
the cloud is shredded and incorporated into in the hot, 
post-shock gas. 

In our study of the effects of metallicity for a given 
shock velocity, we find that only in the slowest shocks 
does an increase in metallicity consistently lead to a de- 
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Fig. 5. — Dust mass evolution (surviving mass) for the nine dust species tracked in this work versus time in units of cloud-crushing time. 
The line styles correspond to the four different shock velocities, 10 3 km s _1 (solid), 3 x 10 3 km s _1 (dot-dashed), 5 x 10 3 km s _1 (dashed), 
and 10 4 km s _1 (dotted). The colors correspond to the three different metallicities, Z = 1 Zq (blue), 10 Zq (red), 100 Zq (green). Panels 
(b), (c), and (h) do not fully span the range to 1 on the y-axis. 



vious work. In the most extreme case, only 24% of the 
Fe dust mass survived. If we look at dust survival for 
two other commonly studied grain species for the same 
extreme case, 44% and 7% of the initial dust mass re- 
mains for carbonaceous (C) and silicate (SiC^) grains, 
respectively. 

For example, consider the differences in dust survival 
across grain species. The difference, ^survival? 

between 

most destroyed grains and least destroyed grains ranges 
from as high as 96% in Simulations 4 and 5 (AI2O3 com- 
pared to Si) to as low as 45% in Simulation 12 (again, 
AI2O3 compared to Si). This variation stems primar- 
ily from the initial distributions in grain radius. A high 
degree of destruction is observed in the species with a 
majority of their mass locked up in small grains. De- 
struction is reduced significantly for those species with 
appreciable mass in large grains. As noted in Paper 1, 
grains with initial radii less than 0.1 /im are easily oblit- 



erated, while larger grains take a considerable duration 
to be substantially sputtered. If the actual size distri- 
butions of dust grains in SNRs were greatly disparate 
from the distributions assumed here, the final dust mass 
survival fractions might change considerably. 

4.2. Density- Temperature Phase-Space 

In an attempt to understand the differences in the 
physical environments between the twelve simulations, 
we present phase-space diagrams of density and temper- 
ature to illustrate where the dust particles spend their 
time (Figure [6]). We also over-plot erosion rate contours 
for C grains as a function of density and temperature 
to determine where in phase-space the majority of the 
dust destruction occurs. The most apparent difference 
between simulations is the rise in gas temperatures with 
increasing shock velocity, a direct result of shock- heating. 
The increase in temperature pushes the gas into a regime 
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TABLE 2 
Final Dust-mass Survival Fraction 



Simulation 


AI2O3 


C 




1 


0.039 


0.947 


0.879 


2 


0.184 


0.953 


0.902 


3 


0.661 


0.964 


0.938 


4 


0.001 


0.714 


0.369 


5 


0.001 


0.714 


0.371 


6 


0.118 


0.659 


0.409 


7 


0.001 


0.610 


0.181 


8 


0.001 


0.597 


0.171 


9 


0.000 


0.531 


0.137 


10 


0.002 


0.585 


0.120 


11 


0.002 


0.568 


0.113 


12 


0.001 


0.442 


0.067 


Simulation 


MgSi0 3 


Si0 2 


MgO 


1 


0.413 


0.910 


0.827 


2 


0.518 


0.931 


0.857 


3 


0.789 


0.951 


0.922 


4 


0.038 


0.411 


0.347 


5 


0.039 


0.413 


0.351 


6 


0.162 


0.441 


0.399 


7 


0.014 


0.203 


0.184 


8 


0.014 


0.192 


0.177 


9 


0.007 


0.156 


0.152 


10 


0.018 


0.128 


0.132 


11 


0.016 


0.121 


0.126 


12 


0.008 


0.072 


0.084 


Simulation 


Si 


FeS 


Fe 


1 


0.998 


0.643 


1.000 


2 


0.998 


0.708 


1.000 


3 


0.993 


0.816 


0.989 


4 


0.964 


0.165 


0.959 


5 


0.960 


0.145 


0.940 


6 


0.815 


0.173 


0.693 


7 


0.849 


0.055 


0.801 


8 


0.822 


0.044 


0.752 


9 


0.660 


0.007 


0.465 


10 


0.726 


0.031 


0.603 


11 


0.688 


0.025 


0.547 


12 


0.446 


0.005 


0.235 


Note. — For simulation parameters, refer to 



Table Q] 



of higher grain erosion, especially when the gas is still 
compressed to relatively high densities. We also note 
that the large number of time steps at high temperature 
and low-density seen across all simulations is a result of 
the simulation end-state, at which point the cloud has 
been shock-heated and shredded. 

In reviewing some of the subtle differences between 
simulations, we identify features that help to explain 
the differences between final dust masses. First, we be- 
lieve the high number of time steps spent by particles 
at high densities but very low temperatures in Simu- 
lation 3 (v s hock = 10 3 km s _1 , Z = 100 Zq) is likely 
the dominating factor in the decreased dust destruction 
when compared to the lower metallicity simulations (1 
and 2). In this scenario, the temperatures of the gas 
drop low enough from the radiative cooling to slow grain 
erosion. Although Simulation 6 (v s hock = 3 x 10 3 km s _1 , 
Z = 100 Zq) shows a similar population of high-density, 
low-temperature dust particle environments, it is not 
sufficient to produce the same consistent decrease in 
dust destruction, owing primarily to the overall increase 
in gas temperature. The pile up of dust particles at 
T = 1000 K, seen in Simulations 3 and 6, is a d irect 
result of the temperature floor described in Section [2731 



For 

^shock — 5 x 10 3 km s x , there is a distinct feature 
present in Simulation 9 with Z = 100 Zq that separates 
it from the lower metallicity cases. As a result of the en- 
hanced metal cooling, the gas is able to cool and condense 
beyond a hydrogen number density of — 10 4 cm -3 , 
but it remains within a temperature regime to allow for 
appreciable grain erosion, explaining the increased dust 
destruction observed in the dust mass evolution. 

In the case of the highest velocity shock, the intense 
heating appears to overwhelm the effects of radiative 
cooling, even for the highest metallicity simulations, as 
the phase plots for Simulations 10-12 appear nearly iden- 
tical in the qualitative sense. Given the similarity be- 
tween the physical environments experienced by the dust 
particles, it would appear that in this regime the domi- 
nant cause for increased dust destruction at high metal- 
licity is the increase in the erosion rates. This is not an 
unreasonable explanation, since the dust particles spend 
the majority of their time above 10 6 K, the point at which 
the dust destruction rates begin to vary significantly for 
higher metal abundances. 

5. SUMMARY AND DISCUSSION 
The above results can be summarized as follows: 

• The degree of destruction varies widely across the 
explored parameter space. For the least destruc- 
tive cases (low shock velocity, high gas metallicity), 
we see near-negligible grain destruction for some 
species (Fe and Si) and the highest survival rates 
for all other species. In the maximally destructive 
cases (high shock velocity, high gas metallicity), 
we find complete to near-complete destruction for 
multiple grain species (AI2O3, FeS, and MgSiOs). 

• The relative velocity between the reverse shock and 
the cloud of eject a material is a dominant factor 
in determining the survival rate of dust grains; an 
increase in velocity correlates directly with an in- 
crease in dust destruction. The most extreme ex- 
ample is for silicate (SiC^) grains, where an ad- 
ditional 88% of the initial dust mass is destroyed 
between the simulations with v s h oc k = 10 3 km s _1 
and those with v s hock = 10 4 km s _1 . 

• There appears to be a threshold shock velocity be- 
tween 1000 and 3000 km s _1 at which the amount 
of dust destruction at early times (t < t cc ) is in- 
creased considerably. 

• In the physically expected, highly metal-enriched 
simulations, we find ranges of dust mass survival 
fraction of 44-96%, 7-95%, and 24-99% for the of- 
ten studied grain species of C, Si02, and Fe, re- 
spectively. 

• As in Paper I, the initial distribution in grain radius 
greatly influences the overall survival of a given 
species, with grains with initial radii below 0.1 /im 
being destroyed in most simulations. 

• For the high- metallicity (Z = 100 Zq) simulations, 
the balance between the increased radiative cooling 
and increased sputtering yield depends strongly on 
the speed of the shock. When the shock velocity is 
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log[n H ] 



Fig. 6. — Phase plots of hydrogen number density (tih in cm 3 ) vs. temperature (T in K). The color scale represents the total number 
of times steps (N) spent by all dust particles in the simulation. The contour lines show the dust destruction rate, \og(da/dt), in units of 
fim yr _1 for C grains. Each panel is one of the 12 simulations, labeled by the number in brackets [Table [l]. 
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below 3000 km s _1 , the increased cooling lowers the 
temperature of the gas enough to suppress dust de- 
struction. At high velocities, the gas temperature 
is driven up to a high-erosion regime. 

In comparing this new suite of simulations to that of 
Paper I, the key differences are the higher metallicity 
cooling functions, the erosion rates for enhanced metal 
abundances for all elements up to iron, the exploration 
of higher shock velocities (v s hock = 10 4 km s _1 ), and 
the suppressed cooling in the ambient medium surround- 
ing the cloud. As a consistency check, we find that 
the simulations with Z = Z@ for shock velocities of 
^shock = 1000 km s~\ 3000 km ~\ and 5000 km s" 1 
produce very similar dust mass survival fractions, indi- 
cating that these changes do not invalidate our previous 
results. In fact, the degree of dust survival does not begin 
to deviate significantly from Paper 1 until the metallicity 
is increased to Z = 100 Z , as expected. 

We also find that our dust survival fractions remain 
in agreement with observations of IR emission in SNRs. 
Specifically, the continued survival of Si grains matche s 
the high abundances found in Cas A by Rho et al.| ([2008), 
while complete destruction of AI2O3 grains and survival 
of amorphous c arbon grains match obse rvations of SNR 



1E0102-7219 by Sandstrom et al. (2009) 



Ideally, we would have liked to explore a metallicity 
regime in which the gas in the eject a cloud is saturated 
with metals, effectively mimicking clouds of pure oxygen, 
silicon, sulfur, argon, or calcium. Using oxygen as an ex- 
ample, an abundance ratio of order one would require a 
metallicity of Z ~ 1000 Z@ since the solar abundance 
ratio is no/nu ^ 5 x 10 -4 . Unfortunately, due to the 
extremely short cooling times, this proved computation- 
ally prohibitive. The simplest solution is to limit the 
simulation time-step by the cooling time and allow the 
simulation to run for a longer physical time, but it is also 
the most computationally expensive. In order to probe 
the highest metallicity regimes, a more resource-efficient 
method would need to be implemented without sacrific- 
ing the validity of the simulation. Such a solution is left 
to future work. 

Currently, we still operate in the limit of strictly ther- 
mal sputtering, since our dust grains are directly coupled 



to the flow of the gas. Because large grains would require 
a greater transfer of momentum from the gas before be- 
coming completely entrained in the flow, it is possible 
that the dominant mode of sputtering for such grains 
would be in the non-thermal limit, as high velocity gas 
washes over them. In order to decouple the dust parti- 
cles from the motions of the gas, we would have to track 
separated grain populations of varying radii. The num- 
ber of radius bins required to accurately track the entire 
dust mass for all nine grain species is uncertain, but it 
could be computationally significant. As mentioned in 
Paper I, it may not be unreasonable for the dust parti- 
cles to be coupled to the flow if the grains are charged and 
magnetic fields permeate the ejecta material. However, 
magnetized ejecta clouds would require an investigation 
of the effects of b etatron acceleration on charged grains 
([S^[l^[l978|. 

Finally, we investigated the possibility of transition- 
ing from post-processing of our dust particle histories to 
compute erosion to an "on-t he-fly" method that would 
sputter the grains as the simulation ran. This would al- 
low metals that are sputtered off the grains to feed back 
into the ambient medium and increase the metallicity. 
However, if we assume that the gas within the ejecta 
cloud is already heavily metal-enriched, the amount of 
metals released into the gas in this way would not al- 
ter the metallicity enough to make a substantial differ- 
ence. As we have shown, it takes a large change in the 
metal abundance to create an significant difference in the 
amount of dust that survives the shock-cloud interaction. 
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